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Abstract 

Dijet events in photon-proton collisions in which there is a large pseudorapidity separa- 
tion Ar] > 2.5 between the two highest By jets are studied with the HI detector at HERA. 
The inclusive dijet cross sections are measured as functions of the longitudinal momentum 
fractions of the proton and photon which participate in the production of the jets, xj, ets and 
xif ets respectively, Arj, the pseudorapidity separation between the two highest Et jets, and 
E^ p , the total summed transverse energy between the jets. Rapidity gap events are de- 
fined as events in which Ej^ p is less than E^, for E^ 1 varied between 0.5 and 2.0 GeV. 
The fraction of dijet events with a rapidity gap is measured differentially in Arj, Xp £ts and 
xiy ets '. An excess of events with rapidity gaps at low values of E™ 1 is observed above the 
expectation from standard photoproduction processes. This excess can be explained by the 
exchange of a strongly interacting colour singlet object between the jets. 
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Figure 1: The diffractive deep inelastic scattering process (a) compared with the rapidity gaps 
between jets process (b) at HERA. 

1 Introduction 



Events with large rapidity gaps in the hadronic final state coupled with the presence of a hard 
scale have long been recognized as an ideal place to study the interplay between long and short 
distance physics [[]]]. Measurements of diffractive deep inelastic scattering at HERA 0,|3|], for 
example, have led to a phenomenology which mixes the languages of perturbative quantum 
chromodynamics (QCD) and Regge theory. In this case the proton remains intact or dissociates 
into a low mass system: the momentum transfer across the gap is small. In other words there 
is no hard scale at the proton vertex (see figure [TJ(a)). A special class of events in which there 
is a large momentum transfer across the gap has been observed in proton - antiproton collisions 
at the Tevatron [0-0] and is the subject of this study. Such events are characterised by two 
high E T jets in the final state separated by a large rapidity interval and little or no energy 
flow between the jets (see figure |l](b)). At such large momentum transfers the Regge-inspired 
phenomenology of diffractive deep inelastic scattering predicts cross sections that are orders 
of magnitude smaller than those measured. In a pp collision, the process may be visualised as 
the hard scattering of one parton from each proton via the exchange of a strongly interacting 
colour singlet object, the "perturbative pomeron". At HERA the ZEUS Collaboration observed 
this process in photoproduction events, in which a quasi real photon emitted from the incoming 
positron interacts with the proton The photon can interact either as a point-like object, 
which in leading order QCD is termed the direct process, or it can fluctuate into a hadronic state 
from which a parton can interact with a parton in the proton, the resolved process. In resolved 
events, parton - parton hard scattering via the exchange of a colour singlet object should be 
present, and it is one of the aims of this paper to search for these events. 

It might be expected that the rate for such a process should be calculable using only pertur- 
bative methods, since the presence of a large momentum transfer ensures that the gap production 
mechanism is dominated by short distance physics QTTp: the perturbative pomeron couples to 
individual partons. The presence of a rapidity gap together with the large rapidity separation 
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of the jets and large momentum transfer means that one is in a perturbatively calculable Regge 
limit of QCD, i.e. s ^> — i ^> A.q CD , where t is the momentum transfer across the gap and 
s is the parton - parton centre of mass energy. Calculations using the leading logarithmic ap- 
proximation (LLA) of BFKL JT^-[T4|] have been performed [|l5l,|l6|]. They have been found to 
describe an observed excess of dijet events with low levels of hadronic activity between the jets, 
over the expectation from standard QCD processes, at the Tevatron [|T3, IT)]. The situation is 
complicated, however, by the possibility that interactions between spectator partons in the col- 
liding hadrons can destroy the gap [|18l-|2"l"|] . This introduces a potentially large non-perturbative 
component into such calculations, making absolute predictions of both rates and differential 
distributions problematic. It has been noted by several authors [[T^,^2[^3]] that this problem can 
be avoided, and perturbative predictions are in principal possible, if a rapidity gap is defined 
in terms of the energy E%? p 3> Aqcd in the pseudorapidity region between the jets. It is the 
aim of this paper to search for dijet events with rapidity gaps in a way which will allow such 
theoretical calculations to be compared to the experimental data in the future. An alternative 
approach to BFKL calculations is offered by the Colour Evaporation Model, in which the gap 
formation mechanism itself is considered to be purely non-perturbative. This model has also 
been shown to describe the Tevatron and HERA data [ |24] ] . 

After a brief description of the HI detector in section [Z[ the definition of a rapidity gap event 
and the selection of the event sample are discussed in section 0. The Monte Carlo models used 
to correct the data and to compare to the BFKL predictions for the rapidity gap production rates 
are described in section In section |5J the data are presented, and the results are discussed in 
section 151. 



2 The HI Detector 

A detailed description of the HI apparatus can be found elsewhere [E5J|. The following briefly 
describes the detector components relevant to this analysis. 

A liquid argon (LAr) calorimeter covers the range in polar angle 4° < 9 < 153° (3.35 > 
T] > —1.43) with full azimuthal coverage^. The LAr calorimeter consists of an electromagnetic 
section with lead absorbers and a hadronic section with steel absorbers, of combined depth 
between 4.5 and 8 interaction lengths. Both sections are highly segmented in the transverse and 
longitudinal directions with about 44000 cells in total. The absolute hadronic energy scale is 
known to 4% for this analysis. The polar region 153° < 9 < 177.8° (-1.43 > r] > -3.95) is 
covered by the SPACAL a lead- scintillating fibre calorimeter with both electromagnetic 
and hadronic sections, with a combined depth of 2 interaction lengths. The hadronic energy 
scale uncertainty is presently 7%. Tracking information is provided by the two concentric drift 
chambers of the central tracker, covering the pseudo-rapidity range —1.5 < r] < 1.5, which lie 
inside a 1.15 T solenoidal field. 

The luminosity is measured from the reaction ep — > ejry with two TICl/TlBr crystal calorime- 
ters p5|], the electron and photon taggers, installed in the HERA tunnel. The electron tagger is 
located at z = —33 m from the interaction point in the direction of the outgoing lepton beam 
and the photon tagger at z = —103 m. 

l 8 is measured relative to the outgoing proton beam direction which defines the positive z axis and the forward 
direction. Pseudo-rapidity is defined as r; = — ln(tan#/2). 
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3 Event Selection and Kinematic Reconstruction 



The data for this analysis were collected with the HI detector during the 1996 running period, 
when HERA collided 27.6 GeV positrons with 820 GeV protons. An integrated luminosity 
of 6.6 ptT 1 is used. Photoproduction events were selected by detecting the scattered positron 
in the electron tagger of the luminosity system. This restricts the virtuality of the photon to 
Q 2 < 0.01 GeV 2 . The photon-proton centre of mass energy was restricted to the range 165 < 
W < 233 GeV to avoid regions of low electron tagger acceptance. 

Two triggers were used in the analysis, one based on tracking requirements in the central 
tracker, the other on energy in the electromagnetic section of the SPACAL calorimeter. The 
combined efficiency of the two triggers, was greater than 60% for all measured data points. 
Both triggers also required an electron to be detected in the electron tagger. 

Hadronic final state objects were defined using a combination of tracking and calorimetric 
information. An algorithm was used which avoids the double counting of tracks and calorimeter 



clusters [J27|]. The inclusive k± clustering algorithm was applied to the hadronic final state 
objects. The algorithm was run in the 'covariant px scheme' with the R parameter set to 1 
[[2~S|,[2~9"|]. In this inclusive mode, the k±_ algorithm merges every final state object in the event, 
excluding the scattered electron, uniquely into a list of massless jets, ordered in E T . Events 
were kept if the two highest E T jets in the event satisfied the following criteria: 

E^' 1 > 6.0 GeV (1) 
E J T et ' 2 > 5.0 GeV (2) 

T fet,l T] jet,2 < 2 65 (3) 

2.5 < A V = \rf et ' 1 - rf et ' 2 \ < 4.0 (4) 
The first two cuts are asymmetric to avoid regions in which NLO QCD predictions become 



unstable [[30]] , in order to facilitate the comparison of the data with such calculations when they 
become available. The third cut ensures that the most forward jet is well contained within the 
LAr calorimeter. The fourth cut ensures a large jet-jet separation. No backward i] cut on the jets 
was necessary, as the other kinematic cuts force the backward jet to be within the acceptance of 
the HI detector. 

The total transverse energy between the two highest E T jets, E^ ap , is defined as the sum of 
the energy of all jets whose pseudorapidities, rf et ' % , lie between their axes, i.e. 

pgap _ \ ^ F jet,i jet jet,i jet 

T — 2.^1 T J 'I forward ^ '/ ^ '/backward' 

i>2 

In addition to the kinematic variables defined above, two variables x^ ts and x p ets are de- 
fined as the fractional longitudinal momentum of the photon and proton participating in the 
production of the two highest E T jets 

lets = *^=1,2 ^ flij (6) 

X ? - 2E p U) 

where E p is the energy of the incoming proton and the sum in the denominator of the xl^ ts cal- 
culation runs over all hadronic final state objects in the event, excluding the scattered electron. 
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Monte Carlo 


a s 


Proton PDF 


7 PDF 


Pf m (GeV) 


P™ p (GeV) 


PYTHIA 5.7 


1-loop 


GRV-LO @ 


GRV-LO [|5|] 


2.2 


1.5 


HERWIG 6.1 + JIMMY 


2-loop 


GRV-LO 


GRV-LO 


1.8 


1.8 



Table 1: Monte Carlo parameters. 

4 Monte Carlo simulations and data corrections 

The Pythia 5.7 [BTJ] and Herwig 6.1 [|32p Monte Carlo event generators were used to correct 



the data for detector acceptance and bin migration effects, and for model comparisons. Both 
generators simulate the direct and resolved production of dijets by quasi real photons. The 
hard scattering matrix elements are calculated to leading order, regulated by a cut-off, P™ in . In 
addition to the primary hard scatter, both generators contain models to simulate the effects of 
multiple parton-parton interactions in a resolved photoproduction event. In Pythia, the prob- 
ability to have several parton-parton collisions in a single event is modeled using the parton 
densities and the usual leading order matrix elements [TJITJ. As for the primary hard scatter, 
the matrix elements are divergent and must be regulated by a cut-off, P™ p , which is the main 
free parameter in the model. This parameter was tuned to give the best description of the HI 
data, after all kinematic cuts, used in this analysis^. In the case of Herwig, multiple interac- 
tions were simulated by the Jimmy package [J53|]. The main free parameter in the Herwig + 
Jimmy model is P™ p , as for Pythia, which in the current version of Jimmy must be set equal 
to P™ m , the cut off for the matrix elements for the primary hard scatter. Again, this parameter 
was tuned to give the best description of the HI data, after all kinematic cuts, used in this anal- 
ysis. The settings of the above parameters and parton densities used in each of the Monte Carlo 
simulations are shown in table []]. 

For both Herwig and Pythia the direct and resolved photoproduction processes were 
generated separately and added together according to their generated cross sections. The overall 
normalisation of the Pythia sample was scaled by a factor of 0.7 in order to fit the measured 
inclusive dijet photoproduction cross sections, shown in figure |]. Similarly, the Herwig sample 
was scaled by a factor of 1.2. 

Herwig incorporates the BFKL LLA colour singlet exchange cross section for the elastic 
scattering of two partons as computed by Mueller and Tang pT3p. In the limit Ar/ ^> 1 the cross 
section for quark - quark scattering may be approximated as^ 

da(qq^qq) 4 2vr 3 exp(2u y) 

— « {C F a s ) (8) 

dt t 2 {7a s C A ({3)yy 

where 

y = A V = \n f (9) 



2 The simplest multiple interaction model in PYTHIA was used, corresponding to setting switch MS TP ( 82 ) =1. 
3 We used a version of HERWIG in which the BFKL LLA colour singlet exchange cross section is valid for all 



A77 0. 
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and 



u; = CU(41n2/7r)a s . (10) 

Here, 1 + cu is the perturbative pomeron intercept, Cp = | is the usual colour factor for quark- 
quark scattering, Ca is the number of colours and ( is the Riemann ^-function. The values of 
a s in equations ^] and [TT] are free parameters in the LLA, and are each chosen to be 0.18. This 
corresponds to a choice of pomeron intercept of 1.48. 

BFKL pomeron exchange has not yet been implemented in Pythia. Colour singlet ex- 
change events were modeled by high-t photon exchange, with a scale factor of 1200 applied to 
the generated cross section. The reason for using this process is twofold; firstly, to obtain a sam- 
ple of events from Pythia in which the hadronic final state is that of a colour singlet exchange 
process, for the purposes of detector acceptance and bin migration corrections. Secondly, the 
high-t photon exchange process provides a useful way of ascertaining how sensitive the data are 
to the underlying dynamics of colour singlet exchange, since the dynamics of photon exchange 
are very different to those of BFKL pomeron exchange. This will be discussed further in section 

The Pythia Monte Carlo sample described above, including high-t photon exchange, was 
passed through a full simulation of the HI detector and used to correct the data for detector 
acceptance and bin migration effects. The scaled high-t photon exchange sample was added 
to the scaled photoproduction sample to fit the measured excess of rapidity gap events. The 
data were also corrected using the Herwig + BFKL sample, and the difference assigned as a 
systematic error. For the purposes of detector corrections, the BFKL sample was scaled by a 
factor of 0.8 and then added to the scaled photoproduction sample, again to fit the measured 
excess. 

The transverse energy flow around the jet axes predicted by the Herwig and Pythia gen- 
erators is compared to the uncorrected dataf] in figure [|, for all events surviving the selection 
cuts described in section |3|, after being passed through a complete simulation of the HI detector. 
Both Pythia and Herwig give a good description of the data both inside and outside the jets. 



5 Results and Model Comparisons 

The ep inclusive dijet cross sections, corrected for detector effects, in the kinematic range de- 
scribed in section [3], are given in table [2] and shown in figure [3[ The cross sections are defined 
at the level of stable hadrons. The inner error bars show the statistical error and the outer error 
bars show the statistical and uncorrected systematic errors added in quadrature. The dominant 
uncorrected systematic errors arise from the uncertainty in the trigger efficiency (~ 5%) and 
the model dependence of the Monte Carlo correction procedure (~ 5%). The correlated sys- 
tematic errors are shown as the solid band below the plots. The dominant error comes from the 
uncertainty in the hadronic energy scale of the LAr calorimeter. In all plots, the bin sizes are 

4 The r\ profile shows Stj = ^ ob j cct — rf ct , weighted by the object transverse energy, for objects within 1 radian 
in <fi of the centre of the jet. The <f> profile is similarly defined for objects within 1 unit of pseudorapidity of the jet 
centre. 
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chosen to keep the effects of migrations to an acceptable level. In figure ^] the Herwig and 
Pythia Monte Carlo generator predictions are shown without additional colour singlet ex- 
change components added. 

Figure 0(a) shows the dijet cross section differential in E^ ap . For E^ ap < 0.5 GeV there 
is a marked excess in the data over the prediction of both Monte Carlo generators although 
Herwig predicts a larger cross section than Pythia. Figure |](b) shows the dijet cross section 
differential in A77. Both Pythia and Herwig give reasonable descriptions of the shape of the 
distribution. In figure [|(c), the dijet cross section differential in x^ ts is plotted. Both generators 
fail to describe the shape of the x^ ts distribution. Figure |3](d) shows the dijet cross section 
differential in x p ets . The minimum kinematically accessible x p ets is set by the requirement that 
Ar] > 2.5, i.e. the parton - parton centre of mass energy must be largeQ. Both Herwig and 
Pythia describe the shape of the x p ets distribution well. 

In order to look in more detail at the excess of events with low values of E^ ap it is helpful 
to look at the gap fraction, primarily because the bulk of the systematic errors cancel. The gap 
fraction is formed by taking the ratio of the dijet cross section, with the additional requirement 
that E 9 T ap < ESp*, to the inclusive dijet cross section, for Ef l = 0.5, 1.0, 1.5 and 2.0 GeV. 
The motivations for choosing several different values of E™ 1 to define the gap fractions are 
primarily theoretical. As discussed in section [l], perturbative calculations of the gap fraction are 
possible if E^ is chosen to be sufficiently large P2],p3[], specifically y— 1 ^> E^ 1 Aq C d- 
Such calculations are able to predict the change in the gap fraction as E™ 1 is increased. The 
cross sections for events with E^ ap < E™ 1 are given differentially in A77, x^ ets and x J p ets in 
tables D, |] and |] respectively, along with the associated gap fractions. 

In figure |](a) the gap fraction differential in Arj is shown for the 4 different choices of 
E™*. Also shown are the predictions of Herwig and Pythia. There is a clear excess over the 
Pythia prediction which persists up to the largest value of E™ 1 and increases with Ar]. One 
would naively expect the gap fraction to fall exponentially with increasing A77 in the absence of 
a colour singlet exchange component, given the assumption that multiplicity fluctuations in the 
hadronic final state obey Poisson statistics. This expectation is borne out by the Pythia Monte 
Carlo. The data do not display this behaviour, and possibly rise in the highest Ar] bin. This 
behaviour is indicative of the presence of a colour singlet exchange component in the data. 
The shape of the Herwig distribution is much flatter than that of Pythia, and does not fall 
exponentially as Ar] increases. Herwig is closer to the data than Pythia, although at large 
Ar] there is again a noticeable excess in the data for all E™ 1 . The difference between the 
Herwig and Pythia predictions is due to the different models of hadronisation in the two 
generators: Jetset [ [3T| ] in the case of Pythia, and the cluster fragmentation model in Her- 
wig [TJ2p. At the level of parton showering, i.e. pre-hadronisation, both generators exhibit the 
expected exponential fall of the gap fraction with increasing Ar]. 

In figure Q(b) the gap fraction for Er^ p < 1.0 GeV is reproduced together with the pre- 
dictions of Herwig and Pythia with colour singlet exchange models added. In this and 
subsequent figures, the BFKL cross section was not scaled to fit the data. The normalisation 
was solely determined by the choice of a s = 0.18 as described in section |j. The Pythia high- 
\t\ photon exchange prediction was scaled by a factor of 1200. It is worth repeating here that 

5 The parton level variables x p and x 1 are correlated, since the rapidity separation between the outgoing patterns 

y = hx(x p XjW 2 / - i). 



9 



the Pythia model is not intended to be a candidate for the observed excess - it is not a strongly 
interacting process. The rationale behind including the curve is to test the sensitivity of the 
data to the underlying dynamics of the exchange. In particular, photons couple only to quarks, 
whereas a gluonic object such as the BFKL pomeron couples preferentially to gluons. As is 
evident from the figure, there is no significant difference in the shapes of the Herwig and 
Pythia distributions: both models are compatible with the data. 

In figure |5|(a) the gap fraction is plotted differentially in x'i^ ts . For E%? p < 0.5 GeV, there 
is an excess visible in the data for x^ ts < 0.75 over the Monte Carlo predictions, although 
as noted above, Herwig predicts a larger cross section at low values of E^ p than Pythia. 
The gap fraction rises significantly at high values of x^ ts . This effect is reproduced by both 
generators, and is due to the fact that, in leading order QCD, direct photoproduction events (high 
x J y ets ) have quark propagators between the outgoing partons associated with the jets, whilst the 
majority of resolved events (low a£ eis ) have gluon propagators. Quark exchange diagrams lead 
to a lower probability of radiation in the rapidity region between the jets than gluon exchange 
diagrams, and hence the gap fraction increases at large x^ ts . 

In figure |5](b) the gap fraction for E^ ap < 1.0 GeV is reproduced together with the predic- 
tions of Herwig and Pythia with colour singlet exchange models added. A better description 
of the gap fraction for x^ ets < 0.75 is achieved for both models, but at the expense of too high 
a gap fraction at large x^ ts . 

In figure ^(a) the gap fraction differential in x p ets is shown. Again, there is an excess over 
the Pythia prediction which persists up to the largest value of E^ 1 . In this case, the excess is 
present in all bins. The tendency for the data and Monte Carlo predictions to rise at low x p ets is 
due to the correlation between x p ets and x^ ets : for large At/, small x p ets must be compensated 
by large x^ ts . The excess over Herwig is less pronounced. In figure |^(b), the gap fraction for 
E^ ap < 1.0 GeV is replotted with the predictions of Herwig and Pythia with colour singlet 
models added. Both models are able to describe the data. 



6 Discussion 

It is clear from the gap fractions presented in figures |J, |5] and |6| that there is a highly significant 
excess of events with low energy flow between the jets over that predicted by the Pythia gen- 
erator before the addition of a colour singlet exchange model. The measurements also show 
an excess over the predictions of Herwig, although the effect is less pronounced. The gap 
fraction in the range 3.5 < At/ < 4.0 for E%? p < 0.5 GeV is approximately 10%. This provides 
the closest comparison with the result of the ZEUS Collaboration, who found a gap fraction of 
approximately 11% at At/ = 3.7, but in a slightly different kinematic range and with a different 
definition of a rapidity gap JT0| ] . 

The large 4-momentum transfer across the rapidity gap forced by the selection of high trans- 
verse energy jets means that standard Regge inspired phenomenology cannot explain this ex- 
cess. One particular solution to this discrepancy has been investigated here, namely the addition 
of a distinct hard colour singlet exchange component to the simulations. The two different mod- 
els of the underlying dynamics of the colour singlet component considered here could not be 
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distinguished by the data. Whilst it is interesting that the LLA BFKL cross section is of the 
right order of magnitude to fit the data for a choice of a s = 0.18, there are significant uncer- 
tainties which make strong conclusions difficult to draw. Equation [8| has been derived using 
the leading logarithmic approximation of BFKL, and the higher order corrections at non-zero 
t are unknown at present. Even within the LLA, there are large uncertainties due to the choice 
of scales, the treatment of the running coupling (which has been treated as a fixed parameter 
in this analysis) and the contribution of higher conformal spin [jTjJ]. A further complication is 
introduced by the issue of rapidity gap survival; colour singlet exchange events will not neces- 
sarily all lead to rapidity gaps. This has not been taken into account in the BFKL predictions 
shown in this paper. 

There is a large uncertainty in the knowledge of the gap formation probability in standard 
photoproduction processes; PYTHIA produces fewer gaps than Herwig + the multiple inter- 
actions package Jimmy. This difference is primarily due to the different hadronisation models 
employed in the two generators. The predictions of PYTHIA, which exhibit an exponential fall 
of the gap fraction as Ar/ increases, are more in line with naive expectations, although there is 
no a priori reason to discount the cluster fragmentation model of Herwig. Finally, even within 
a particular model, the gap formation probability is of course dependent upon the treatment of 
multiple interactions, although the requirement that the level of hadronic activity in the event 
matches the data, as shown in figure [|, provides some contraints. 



7 Summary 

Dijet events in photoproduction have been studied in which there is a large rapidity separation 
between the two highest E T jets. The inclusive dijet cross sections have been measured as 
functions of the longitudinal momentum fractions of the photon and proton which participate 
in the production of the jets, x^ ts and x J p ets respectively, Ar], the pseudorapidity separation 
between the two highest E T jets, and E^ ap , the total summed transverse energy between the 
jets. A significant excess of dijet events with small E^ ap is observed for Ar/ > 2.5 over that 
predicted by standard photoproduction models. 

In order to investigate the excess, the dijet cross sections have been measured with the 
additional constraint that the transverse energy between the jets be less than 0.5, 1.0, 1.5 and 2.0 
GeV. The ratios of these cross sections to the inclusive dijet cross sections, the gap fractions, are 
found to be reasonably well described with the addition of a colour singlet exchange component 
in the form of the LLA BFKL pomeron, although there is little sensitivity in the data to the 
underlying dynamics of the model. 
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1.29 


0.03 


0.19 
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t 0.05 


0.07 


0.07 


X p 


da/dx p ets 
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^uncor 
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0.02 - 0.04 


22.3 


0.5 
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3.9 


0.04 - 0.06 


23.2 


0.6 


2.2 


3.9 


0.06 - 0.08 


08.6 


0.3 
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1.7 


0.08-0.10 
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da/dArj 
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r 
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0.21 
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0.67 


0.02 


0.08 


0.13 
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0.17 


0.01 
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E 9 T ap 
(GeV) 


da/dE^ ap 
(nb/GeV) 
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0.0 - 0.50 
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0.016 
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0.5-1.50 
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3.5-7.00 


0.124 


0.003 


0.014 


0.018 


7.0-12.0 


0.054 


0.001 


0.009 


0.012 



Table 2: The inclusive dijet cross sections differential in xi^ ts , x p ets , An and E^ ap . Also shown 
are the statistical error, 5 stat , uncorrected systematic error, 5 uncor , and correlated systematic 
crror^ ^ cow * 
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0.019 


0.5 


3.1 - 3.5 
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0.006 


0.003 
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0.025 


0.009 


0.010 


0.5 
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0.017 


0.007 
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0.002 
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1.0 
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1.0 


3.5-4.0 
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0.006 


0.004 
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0.029 


0.035 


1.5 


2.5-2.8 
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0.019 
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0.137 


0.010 


0.027 


1.5 


2.8-3.1 


0.157 
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0.020 
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0.014 


0.038 


1.5 


3.1-3.5 
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0.011 
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0.007 
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0.024 
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2.0 
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0.030 


0.178 
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0.033 


2.0 


2.8-3.1 
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0.014 
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2.0 


3.1-3.5 
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0.012 
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2.0 
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0.033 


0.007 


0.006 
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0.034 


0.044 



Table 3: The dijet cross sections differential in A77, with the additional requirement that E^ p < 
Ej^ 1 , shown with the statistical error, S sta t, uncorrelated systematic error, 5 uncor , and correlated 
systematic error, 5 corr . Also shown are the gap fractions, /(A77), defined as the fraction of all 
dijet events with E 9 ^ v < E™*, and their associated statistical and systematic errors. 
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Table 4: The dijet cross sections differential in x^ ts , with the additional requirement that E^ p < 
Etp 1 *, shown with the statistical error, S stat , uncorrelated systematic error, S uncor , and correlated 
systematic error, 5 corr . Also shown are the gap fractions, /(x^ e ' s ), defined as the fraction of all 
dijet events with E^ ap < E^, and their associated statistical and systematic errors. 
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Table 5: The dijet cross sections differential in x J p ets , with the additional requirement that E^ ap < 
E™ f , shown with the statistical error, S stat , uncorrelated systematic error, S uncor , and correlated 
systematic error, 5 corr . Also shown are the gap fractions, f(x p ets ), defined as the fraction of all 
dijet events with E^ ap < E™ 1 , and their associated statistical and systematic errors. 
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Figure 2: Jet profiles represented by the energy flow per jet per unit rapidity (Srj) and per unit 
azimuth (5(f)) around the jet axis for jets forward and backward relative to the proton direction. 
The HI data, uncorrected for detector effects, are shown as points. The solid histogram is the 
prediction of the Pythia simulation, and the dashed histogram that of Herwig, after being 
passed through a full simulation of the HI detector. 
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Figure 3: In (a), the solid points show the dijet cross section differential in E^ p , the summed 
transverse energy between the two highest E T jets, in the kinematic range defined in section 
The inner error bars represent the statistical error, and the outer error bars represent the 
statistical and non-correlated systematic errors added in quadrature. The solid band below the 
plot shows the correlated systematic errors, as described in the text. The dashed line shows the 
prediction of Herwig, scaled by a factor of 1.2, and the solid line that of Pythia, scaled by 
a factor of 0.7. In (b), (c) and (d) the dijet cross sections differential in A77, x^ ts and x J p ets are 
shown. 
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Figure 4: The solid points show the gap fraction differential in A77. The inner error bars repre- 
sent the statistical error, and the outer error bars represent the statistical and all systematic errors 
added in quadrature. Gap events are defined for 4 values of E^ ap , shown in the figures. In (a), 
the gap fractions are compared to the prediction of Herwig (dashed line) and PYTHIA (solid 
line). In (b), the gap fraction is shown for E%? p < 1.0 GeV, and compared to the Herwig and 
PYTHIA predictions with 2 different models of colour singlet exchange added. The dashed line 
shows Herwig + BFKL and the solid line shows Pythia + high-t photon exchange. The 
photon exchange cross section is scaled by a factor of 1200 (see text). 
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Figure 5: The solid points show the gap fraction differential in x^ ts , for Ar) > 2.5. The inner 
error bars represent the statistical error, and the outer error bars represent the statistical and all 
systematic errors added in quadrature. Gap events are defined for 4 values of E%? p , shown in the 
figures. In (a), the gap fractions are compared to the prediction of Herwig (dashed line) and 
PYTHIA (solid line). In (b), the gap fraction is shown for E^ p < 1.0 GeV, and compared to the 
Herwig and Pythia predictions with 2 different models of colour singlet exchange added. 
The dashed line shows Herwig + BFKL and the solid line shows Pythia + high-t photon 
exchange. The photon exchange cross section is scaled by a factor of 1200 (see text). 
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Figure 6: The solid points show the gap fraction differential in x J p ets , for A77 > 2.5. The inner 
error bars represent the statistical error, and the outer error bars represent the statistical and all 
systematic errors added in quadrature. Gap events are defined for 4 values of E^ ap , shown in the 
figures. In (a), the gap fractions are compared to the prediction of Herwig (dashed line) and 
PYTHIA (solid line). In (b), the gap fraction is shown for E^ ap < 1.0 GeV, and compared to the 
Herwig and Pythia predictions with 2 different models of colour singlet exchange added. 
The dashed line shows Herwig + BFKL and the solid line shows Pythia + high-t photon 
exchange. The photon exchange cross section is scaled by a factor of 1200 (see text). 
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